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The coalescence of massive black holes generates gravitational waves (GWs) that will be measur- 
able by space-based detectors such as LISA to large redshifts. The spins of a binary's black holes 
have an important impact on its waveform. Specifically, geodetic and gravitomagnetic effects cause 
the spins to precess; this precession then modulates the waveform, adding periodic structure which 
encodes useful information about the binary's members. Following pioneering work by Vecchio, we 
examine the impact upon GW measurements of including these precession-induced modulations in 
the waveform model. We find that the additional periodicity due to spin precession breaks degenera- 
cies among certain parameters, greatly improving the accuracy with which they may be measured. 
In particular, mass measurements are improved tremendously, by one to several orders of magnitude. 
Localization of the source on the sky is also improved, though not as much — low redshift systems 
can be localized to an ellipse which is roughly a few x 10 arcminutes in the long direction and a 
factor of 2 — 4 smaller in the short direction. Though not a drastic improvement relative to analyses 
which neglect spin precession, even modest gains in source localization will greatly facilitate searches 
for electromagnetic counterparts to GW events. Determination of distance to the source is likewise 
improved: We find that relative error in measured luminosity distance is commonly ~ 0.2% — 0.7% 
at z ~ 1. Finally, with the inclusion of precession, we find that the magnitude of the spins them- 
selves can typically be determined for low redshift systems with an accuracy of about 0.1% — 10%, 
depending on the spin value, allowing accurate surveys of mass and spin evolution over cosmic time. 

PACS numbers: 04.80.Nn, 04.25.Nx, 04.30.Db 



I. INTRODUCTION 
A. Background to this analysis 

Observations have now demonstrated that massive 
black holes are ubiquitous in the local universe. It ap- 
pears that all galaxies with central bulges contain black 
holes whose masses are strongly correlated with the prop- 
erties of the bulge PHH . Hierarchical structure formation 
teaches us that these galaxies assembled over cosmic his- 
tory through the repeated coalescence of the dark mat- 
ter halos in which they reside [3|. Taken together, these 
suggest that coalescences of massive black holes should 
be relatively frequent events, especially at high redshift 
when halo coalescences were common [4] . 

Massive black hole coalescences are extremely strong 
gravitational wave (GW) sources. In the relevant mass 
band — thousands to millions of solar masses — the 
GWs these binaries generate are at low frequency (/ ~ 
10 5 — 1CP 1 Hz) where ground-based GW antennae 
have poor sensitivity due to geophysical and other terres- 
trial noise sources. Measuring GWs from massive black 
holes requires going into the quiet environment of space. 
LISA, the Laser Interferometer Space Antenna, is being 
designed as a joint NASA-ESA mission to measure GWs 
in this frequency band; cosmological massive black hole 
coalescences are among its highest priority targets. By 
measuring these GWs, one can infer the properties of the 
source that generated the waves. Some particularly im- 
portant and interesting properties are the masses of the 
binary's members, their spins, the binary's location on 



the sky, and its distance from the solar system barycen- 
ter. Measuring a population of coalescence events could 
thus provide a wealth of data on the cosmological distri- 
bution and evolution of black hole masses and spins. 

Most analyses of how well binary black hole param- 
eters can be determined by LISA measurements have 
ignored the impact of spin- induced precession 0, S 0] • 
Under such an assumption, subsets of parameters can be 
highly correlated with each other, increasing the errors 
in parameter estimation. One such subset comprises the 
binary's "chirp mass" M, its reduced mass /i, and the 
spin parameters and a (which are written out explic- 
itly in Sec. Ill Ap . These four parameters influence the 
GW phase <E>. As discussed in [1,11, the correlation co- 
efficient between /i and (3 is nearly 1. It is thus difficult 
to "detangle" these parameters from one another in a 
measurement. 

Another such subset consists of a binary's sky position, 
orientation, and distance. To see why these parameters 
are strongly correlated, consider the form of the two po- 
larizations of the strongest quadrupole harmonic of the 
gravitational waveform: 

M 5/3r f )2/3 

h+(t)=2 X (l + cos 2 0cos$(i), (1.1) 

h x (t)=-4- cosisin$(t). (1.2) 

(We work in units with G = 1 = c; a convenient con- 
version factor in this system is 10 6 M@ — 4.92 seconds.) 
The quantity i is the binary's inclination relative to the 
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line of sight: cost = L • n, where L, the direction of the 
binary's orbital angular momentum, defines its orienta- 
tion and n is the direction from observer to source. The 
quantity Dl is the luminosity distance to the source, and 
/(<) = (l/27r)d$/di. 

One does not measure the polarizations h+ and h x di- 
rectly; rather, one measures a sum /im(^) m which the 
two polarizations are weighted by antenna response func- 
tions as follows: 

h M {t) = F + (6 N ,(f) N ,ip N )h + (t) + F x (6 N ,4> N ,ip N )h x (t) . 

(1.3) 

(This equation should be taken as schematic; see III CI for 
a more detailed and definitive description.) The angles 
On and <pN denote the location of the source on the sky 
in some appropriate coordinate system. The angle ijiN, 
known as the "polarization angle," fixes the orientation 
of the component of L perpendicular to the line of sight. 
(In other words, L is fixed by i and ipN-) 

Measuring the phase determines chirp mass with high 
accuracy; the fractional error in M. is often ~ 10~ 3 — 
1CP 4 . As far as amplitude is concerned, the chirp mass 
can be regarded as measured exactly. What remains is to 
determine, from the measured amplitude and the known 
M, the angles 8 N , <f) N , ip N , l, and the distance D L . 

As Eqs. (fTTTj) . (TO)) , and (fO]) illustrate, these five pa- 
rameters are strongly correlated. The motion of LISA 
around the sun 1 breaks these degeneracies to some ex- 
tent — the angles On and <f>N appearing in Eq. (|1.3p can 
be regarded as best defined in a coordinate system tied 
to LISA. As the antenna orbits the sun, these angles be- 
come effectively time dependent. The one- year periodic- 
ity imposed by this motion makes it possible to detanglc 
these parameters. Analyses typically find that the posi- 
tion of a merger event at z ~ 1 can be determined, on 
average, to an ellipse which is 1.5 — 2 degrees across in 
the long direction and 1.5 — 2 times smaller in the short 
direction 2 [f| @, G3 ■ The distance to such a binary can 
be determined to 1% — 2% accuracy on average (less in 
some exceptional cases) [f| 0, E3] • 

B. Black hole spin and spin precession 

The preceding discussion ignores an important piece of 
relativistic physics: the precession of each binary mem- 
ber's spin vector due to its interaction with the spacetime 
in which it moves. In general relativity, the spacetime of 
an isolated object can be regarded as having an "elec- 
tric piece," arising from the object's mass and mass dis- 
tribution, and a "magnetic piece," arising from the ob- 



ject's mass currents and their distribution 3 . Spin preces- 
sion consists of a geodetic term, arising from the parallel 
transport of the spin vector in the gravitoelectric field of 
the other hole, and Lense-Thirring terms, caused by the 
gravitomagnetic field of the other hole. The basic physics 
of gravitomagnetic precession can be simply understood 
by analogy with a similar (and closely related) electro- 
magnetic phenomenon — the precession of a magnetic 
dipole fi immersed in an external magnetic field B. An 
object's spin angular momentum S can be regarded as 
a gravitational "magnetic dipole." When immersed in a 
"gravitomagnetic field," one finds that S feels a torque, 
just as a magnetic dipole fi experiences a torque when 
immersed in magnetic field B. In a binary black hole 
system, the gravitomagnetic field arises from the binary's 
orbital motion and the spins of its members. Precession 
thus includes both spin-orbit (geodetic and orbital grav- 
itomagnetic) and spin-spin effects [12j ]. (The major goal 
of the "Gravity Probe B" experiment is to measure the 
effects of geodetic and spin-spin Lense-Thirring preces- 
sion upon a gyroscope in low Earth orbit 

As the spins precess, they do so in such a way that 
the total angular momentum J = L + S1 + S2 is held 
constant; the orbital angular momentum L precesses to 
compensate for changes in Si and S2- As a consequence, 
the inclination angle 1 and polarization angle ipN be- 
come time varying (as do certain other quantities ap- 
pearing in the GW phase function $). Figure [T] shows 
the so-called "polarization amplitude," defined in Sec. 
Ill CI of the waveform measured by a particular detec- 
tor. Without precession, this quantity is modulated by 
the orbital motion of LISA, helping to provide some in- 
formation about the binary's sky position. The polar- 
ization amplitude also depends on the angles l and ipN, 
so it undergoes additional modulation when precession is 
included. Such precession-imposed time variations quite 
thoroughly break many of the degeneracies which have 
been found to limit parameter measurement accuracy in 
earlier analyses. 

It is without a doubt that black holes in nature spin. 
Observations are not yet precise enough to indicate the 
value of typical black hole spins; the evidence to date 
does, however, seem to indicate that fairly rapid rotation 
is common. For example, the existence of jets from active 
systems seems to require non- negligible black hole spin — 
jets appear to be "launched" by the shearing of magnetic 
field lines (supported by the highly conductive, ionized 
material accreting onto the black hole) by the differential 
rotation of spacetime around a rotating black hole [3, 
0j|. Also, observations of highly distorted iron K-a lines 
— a very sharp flourescence feature in the rest frame of 



1 LISA is being designed as a constellation of three spacecraft 
whose centroid orbits the sun with a period of one year; see 
http://lisa.nasa.gov for further details. 

2 It is worth bearing in mind that the full moon subtends an angle 
of about 30 arcminutes. 



3 This analogy is most apt in the weak field. In that limit, one 
can recast the Einstein field equations of general relativity into 
a form quite similar to Maxwell's equations; see for detailed 
discussion. Though the analogy does not fit quite so well in 
strong field regions, it remains accurate enough to be useful. 
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the emitting iron ions — indicate that this emission is 
coming from very deep within a gravitational potential 
(at radii less than the Schwarzschild innermost radius 
6M) and is smeared by near luminal relativistic speeds 
to boot [ijj]. Though perhaps influenced somewhat by 
selection effects 4 , these pieces of evidence are strong hints 
that the black holes which will form the binaries we hope 
to measure will be strongly influenced by spin. 

The only limit in which spin precession can be ne- 
glected is that in which the spins of the binary's members 
are exactly parallel (or antiparallel) to one another and 
to the orbital angular momentum L. Since the target 
binaries of this analysis are created by galactic merger 
processes, their members will almost certainly have no 
preferred alignment — random spin and orbit orientation 
is expected to be the rule. (This expectation is borne out 
by work [l8| showing that jets in active galaxies are ori- 
ented randomly with respect to the disks of their host 
galaxies.) Taking into account spin precession is thus of 
paramount importance for GW observations of merging 
black hole systems. 

A great deal of work has gone into developing fami- 
lies of model waveforms ( "templates" ) sufficiently robust 
to detect GWs from spinning and precessing binaries, at 
least in the context of measurements by ground-based 
detectors 0, M, HI HI M, HI SI . The key issue in 
this case is that the various modulations on the wave- 
form imposed by the binary's precession smear its power 
over a wider spectral range, making it much more diffi- 
cult to detect at the (relatively) low signal-to-noise ratios 
(SNRs) expected for ground-based observations. Not as 
much work has gone into the complementary problem of 
measuring these waves — examining the impact preces- 
sion has upon the precision with which binary proper- 
ties may be inferred from the waves. To date, the most 
complete and important analysis of this type is that of 
Vecchio (26|. Vecchio focuses (for simplicity) on equal 
mass binaries and only includes the leading "spin-orbit" 
precession term. This limit is particularly nice as a first 
analysis of this problem, since it can be treated (largely) 
analytically (cf. discussion in Sec. Ill B of Ref. [2q). 

Vecchio's work largely confirms the intuitive expecta- 
tion discussed above — the precision with which masses 
are measured is substantially improved; in particular, the 
reduced mass of the system can be measured with several 
orders of magnitude more accuracy. Parameters such as 
the sky location of the binary and the luminosity distance 
are also measured more accurately, but only by a factor 
of 2 - 10. 



4 The systems for which we have constraints on spin are systems 
which are actively accreting, and, ipso facto, those which are 
most likely to be rapidly spinning 11711 . 
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FIG. 1: These figures depict the "polarization amplitude" 
A po i(t) of the signal measured in detector I as a function of 
time. The curves are as follows: solid line, xi — X2 = 0; 
dashed line (nearly overlapping the solid line), xi — X2 = 0.1; 
dotted line, xi = X2 = 0.5; and dot-dashed line, xi — X2 = 
0.9. (x = S/m 2 is the dimensionless spin parameter.) The 
top figure covers the last two years of inspiral. The spinless 
curve has periodicity of one year, corresponding to the motion 
of LISA around the sun. Notice that as spin is introduced, 
the curves become more strongly modulated, with the num- 
ber of additional oscillations growing as the spin is increased. 
By tracking these spin-precession-induced modulations, it be- 
comes possible to better measure parameters like mass and 
sky position and measure spin for the first time. The bottom 
figure shows a close-up of the final months of inspiral. 



C. This analysis 

Our goal here is to update Vecchio's pioneering analy- 
sis by taking the precession equations and the wave phase 
to the next higher order and by performing a broader pa- 
rameter survey (including the impact of mass ratio). By 
taking the precession equations to higher order, we in- 
clude "spin-spin" effects — precessional effects due to one 
black hole's spin interacting with gravitomagnetic fields 
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from the other hole's spin. By taking the wave phase 
to higher order, we include, among other terms, a time- 
dependent spin-spin interaction. Finally, when the mass 
ratio differs from 1, the geodetic spin-orbit term causes 
the two spins to precess at different rates, even without 
the spin-spin corrections. 

Including these effects means that we cannot model the 
precession with a simple, analytic rule — we are forced to 
integrate the equations of precession numerically as in- 
spiral proceeds, incurring a significant performance cost. 
Fortunately, the basic "engine" on which this code is 
based [6| runs extremely fast, thanks largely to the use 
of spectral integrators (which, in turn, is thanks to a 
suggestion by E. Berti [7j), so total run time remains 
reasonable. 

The cost in efficiency due to the inclusion of higher- 
order effects is offset by the more complete description 
of the signal they provide. An important consequence is 
that it now becomes possible from GW measurements to 
determine the spin of each member of the binary. With 
Vecchio's approximations, only three components of the 
black holes' vector spins can be determined — enough to 
constrain, but not determine, their spin magnitudes. Our 
more general approach allows us to measure all six vector 
spin components. To our knowledge, this is the first anal- 
ysis indicating how well spin can be measured from merg- 
ing comparable mass binary systems. (As Barack and 
Cutler have shown [27|, spin is very well determined by 
measurements of GWs from extreme mass ratio binaries 
— those in which the system's mass ratio m^/mi < 10~ 4 
or so.) 

Our error estimates are computed using the maximum 
likelihood formalism first introduced in the context of 
GW measurements by Finn 28]. A potential worry is 
that we are using a Gaussian approximation to the like- 
lihood function. This approximation is very convenient 
since it allows us to directly compute a Fisher informa- 
tion matrix. Its inverse is the covariance matrix, which 
directly encodes the estimated 1-a errors in measured 
parameters, as well as correlations among different pa- 
rameters. The Gaussian approximation is known to be 
accurate when the SNR is "high enough" [l,[28j]. 

Unfortunately, it is not particularly obvious what 
"high enough" really means. In our case, we are estimat- 
ing measurement errors on 15 parameters 5 — a rather 
fearsome number to fit. The Gaussian approximation al- 
most certainly underestimates measurement error, since 
it assumes the likelihood function is completely deter- 
mined by its curvature in the vicinity of a maximum, 
missing the possibility of a long tail to large error. We 
thus fear that our estimates are likely to be optimistic, 



5 Two masses; 2 angles specifying the initial orientation of the 
binary's orbit; 4 angles specifying the initial orientation of the 
spins; 2 spin magnitudes; the time at which coalescence occurs; 
the phase at coalescence; 2 angles specifying the binary's position 
on the sky; and the distance to the binary. 



especially for events with relatively small SNR. It would 
be quite salubrious to "spot check" a few cases by di- 
rectly computing the likelihood function in a few impor- 
tant corners of parameter space and comparing to the 
Gaussian predictions. This would both quantify the de- 
gree to which our calculations are too optimistic and help 
to determine how large SNR must be for this approxima- 
tion to be reliable. 

In addition to concerns about the Gaussian approxi- 
mation, it must be noted that the waveform family we 
use for our analysis is somewhat limited. We use a post- 
Newtonian description of the GWs from these binaries. 
Since our analysis requires us to follow these binaries 
deep into the strong field where the usual post-Newtonian 
expansion is likely to be somewhat unreliable, it is likely 
that we are introducing some systematic error. In partic- 
ular, the equations of spin precession that we are using 
are only given to the leading order needed to see spin- 
orbit and spin-spin precession effects [291 ]. Higher spin- 
orbit corrections to the equations of motion and preces- 
sion have recently been derived [30jl. as have their im- 
pact on the the waves' phasing [311 ]. Another analysis 
[321 ] has worked out higher-order spin-spin corrections to 
the post-Newtonian metric, from which it would not be 
too difficult to work out equations of motion and preces- 
sion and then the modification to the waves' phase. It 
would be interesting to see what effect the higher-order 
corrections have on these results. 

Finally, it should be noted that the frequency do- 
main expression of the signal which we are using is de- 
rived formally using a "stationary phase" approximation. 
This approximation is based on the idea that the bi- 
nary's orbital frequency is changing "slowly." The orbital 
frequency is thus well-defined over "short" time scales. 
Quantitatively, this amounts to a requirement that the 
time scale on which radiation reaction changes the or- 
bital frequency, T; nsp , be much longer than an orbital 
period, T or b. Precession introduces a new time scale, 
T prcc , the time it takes for the angular momentum vec- 
tors to significantly change their orientations. For the 
stationary phase approximation to be accurate, we must 
in addition require T proc S> Torb, a- somewhat more strin- 
gent requirement than Ti nsp S> Torb- No doubt, a certain 
amount of error is introduced due to the breakdown of 
this condition late in the inspiral. 

Thus, the results which we present here should be 
taken as indicative of how well LISA is likely to be able 
to measure the parameters of massive black hole binaries, 
but cannot be considered definitive. We are confident 
however that the improvement in measurement accuracy 
obtained by taking spin precession into account is robust. 
Specifically, we see that errors in masses are reduced dra- 
matically, from one to several orders of magnitude. Er- 
rors in sky position and distance are also reduced, but 
by a smaller factor. Such improvement may nonetheless 
critically improve the ability of LISA to interface with 
electromagnetic observatories [U H3, HH]. Finally, the 
added information in the precession signal allows us to 
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measure the spins of the holes. These improvements due 
to precession will certainly survive and play an impor- 
tant role even in an analysis which addresses the caveats 
we list above. 



D. Organization of this paper 

The remainder of the paper is organized as follows. In 
Sec. m we discuss the gravitational waveform generated 
by binary black hole coalescence, focusing upon the slow, 
adiabatic inspiral. Section III Al describes the "intrinsic" 
waveform produced by the motion of the orbiting black 
holes as given in the "restricted post-Newtonian expan- 
sion" of general relativity. Section III Bl then describes 
the post-Newtonian precession equations which we use to 
model the evolution of the spins of a binary's members, 
as well as how those precessions influence the waveform. 
Finally, in Sec. Ill Cl we describe "extrinsic" effects which 
enter the measured waveform through its measurement 
by the LISA constellation. 

In Sec. IIIII we summarize our parameter estimation for- 
malism; this section will be largely review to readers fa- 
miliar with the literature on GW measurements. Section 
IIII Al first summarizes the maximum likelihood formalism 
we use to estimate measurement errors. In Sec. IIIIBI we 
then describe the up-to-date model for the noise which 
we expect to accompany LISA measurements. 

Section [TV] presents our results. After describing some 
critical procedural issues in the setup of our calculations 
in Sec. IIV Al we summarize our results for parameters 
intrinsic to the binary (particularly masses and spins) in 
Sec. IIVBI and for extrinsic parameters (particularly sky 
position and luminosity distance) in Sec. IIV CI In both 
cases, we compare, when appropriate, to results from a 
code which does not incorporate spin-precession physics. 
(This code was originally developed for the analysis pre- 
sented in Ref. 0.) The general rule of thumb wc find is 
that the accuracy with which masses can be determined 
is improved by about one to several orders of magnitude 
when precession physics is taken into account. In addi- 
tion, we find that for low redshift (z ~ 1) binaries LISA 
should be able to determine the spins of the constituent 
black holes with a relative precision of 0.1% — 10%, de- 
pending (rather strongly) on the spin value. Likewise, 
wc find improvement in the measurement accuracy of ex- 
trinsic parameters, though not quite as striking — half 
an order of magnitude improvement in source localiza- 
tion and distance determination is a good, rough rule of 
thumb. 

An important consequence of these improvements is 
that LISA should be able to localize low redshift bina- 
ries — using GW measurements alone — to an elliptical 
"pixel" that is perhaps a few x 10 arcminutes across in 
its widest direction and a factor of 2 — 4 smaller along its 
minor axis. For higher redshift binaries (z ~ 3 — 5), this 
pixel is several times larger, perhaps a few degrees in the 
long direction and tens of arcminutes to a degree or two 



in the narrow one. These results suggest that it should 
not be too arduous a task to search for electromagnetic 
counterparts to a merging binary black hole's GW signal 
[33l [3~il [3^ — particularly at low redshift, these pixel 
sizes are comparable to the field of view of planned large 
scale surveys. 

A concluding and summarizing discussion is given in 
Sec. [V] Along with summarizing our major results and 
findings, we discuss future work which could allow us 
to quantitatively assess the consequences of some of the 
simplifying assumptions we have made. 

At several points in this analysis, we need to convert 
between a source's redshift z and luminosity distance D^. 
To make this conversion, we assume a flat cosmology 
(^totai = 1) with contributions from matter (Qm = 0.25) 
and from a cosmological constant (equation of state pa- 
rameter w = — 1, £1a = 0.75). We also choose a Hubble 
constant Hq = 75 km s _1 Mpc -1 . These choices are in 
concordance with the latest fits presented by the WMAP 
team in their three- year analysis of the cosmic microwave 
background [36( ] . The luminosity distance as a function 
of redshift is then given by 

D (z) = t 1 + ^ c f dz ' (1 4) 

#o Jo y/Sl M (l + z>)3 + n A ' 



II. GRAVITATIONAL WAVES FROM BINARY 
BLACK HOLE INSPIRAL 

The coalescence of a black hole binary can be divided 
into three stages: (1) an adiabatic inspiral, (2) a merger, 
and (3) a ringdown, when the resulting black hole set- 
tles down to its final state. In this paper, we will focus 
on the inspiral. Ringdown waves have been analyzed in 
other work [|| [28|, [33, (38|; the most comprehensive re- 
cent analysis was performed by Berti, Cardoso, and Will 
[3^ ]. The merger waveform, describing the strong field 
and (potentially) violent process of the two black holes 
merging into a single body, has historically been poorly 
understood. Recent breakthr oug hs in numerical relativ- 
ity may soon correct this [40, [4j], [42j • 

The inspiral waveform which will be measured by LISA 
is a combination of the intrinsic waveform created by the 
source and extrinsic features related to its location on the 
sky and modulation effects caused by the motion of the 
detector. In this section we review the relevant physics 
involved in the construction of the waveform. 

For sources at cosmological distances, all time scales 
redshift by a factor 1+z. In the G = c = 1 units that we 
use, all factors of mass enter as time scales; thus, masses 
are redshifted by this 1 + z factor. [Likewise, quanti- 
ties such as spin which have dimension (time) 2 acquire 
a factor (1 + z) 2 , etc.] In the equations written below, 
we do not explicitly write out these redshift factors; they 
should be taken to be implicit in all our equations. When 
discussing results, we will always quote masses as they 
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would be measured in the rest frame of the source, with 
redshift given separately. 



A. Intrinsic waveform 

We treat the members of our binary as moving on qua- 
sicircular orbits; eccentricity is very rapidly bled away by 
gravitational radiation reaction [43[ , so it is expected that 
these binaries will have essentially zero eccentricity by the 
time they enter LISA's frequency band (at least at the 
mass ratios we consider in this paper, 1 < mi/7712 < 10). 
We use the post-Newtonian formalism, an expansion in 
internal gravitational potential U and internal source ve- 
locity v, to build our waveforms. A detailed review of the 
post-Newtonian formalism can be found in the article by 
Blanchet [44J; the key pieces which we will use can be 
found in Refs. Ef El, Ezl 13 ■ 

The post- Newtonian equations of motion, taken to sec- 
ond post-Newtonian (2PN) order, yield the following gen- 
eralization of Kepler's third law relating orbital angular 
frequency f2 and orbital radius (in harmonic coordinates) 



n 2 = ^ 
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(2.1) 



Here M = mi + mi is the total mass of the system, and 
rj = (jl/M, where fi = mim^jM is the reduced mass. L 
is the direction of the orbital angular momentum, and Si 
is the spin angular momentum of black hole i. The mag- 
nitude of the spin can be expressed as Si — Xi m h where 
< Xi < 1- The leading term is the standard result 
from Newtonian gravity. The 0(M/r) term is the first 
post-Newtonian correction; this is the same physics that, 
in solar system dynamics, causes the precession of the 
perihelion of Mercury. The (3((M/r) 3 / 2 ) term contains 
spin-orbit corrections to the equation of motion. Finally, 
the 0((M/r) 2 ) term is a 2PN correction, which also in- 
cludes spin-spin terms. From the equations of motion, 
the orbital energy of the binary E can also be computed; 
see 0. 

The waveform that we will use is the "restricted" 2PN 
waveform. This approximation can be understood by 
writing the waveform (somewhat schematically) as |8| 



4>orb(t) 



(2.2) 



where x labels PN order, m is a harmonic index, and 
^o-[\>{t) = / fl(t')dt' is orbital phase. In the restricted 
post-Newtonian waveform, we throw out all amplitude 
terms except h\ (the "Newtonian quadrupole" term) but 
compute $orb(£) to some specified PN order. The re- 
stricted PN approximation is motivated by the fact that 
matched filtering — matching a signal in noisy data by 
cross-correlating with a theoretical template — is much 
more sensitive to phase information than to the ampli- 
tude. Since the h\ harmonic contributes most strongly to 
the waveform over most of the inspiral, the restricted PN 
approximation is expected to capture the most impor- 
tant portion of the inspiral waveform. It is worth noting, 
however, that there is additional information encoded 
by the harmonics that we are neglecting. Especially for 
the SNRs expected for typical LISA binary black hole 
measurements, this additional information could play an 
important role in measuring source characteristics, as 
pointed out by Hellings and Moore (49j . 

At any rate, within the restricted PN approximation, 
the waveform can be written 



hij (t, x) = - 



4M 5 / 3 {irf) 2 / 3 



C0S$(7j) 

sin <!>(£) 




sin$(t) 
-cos$(t) 




(2.3) 

where |x| is the distance to the binary 
is the "chirp mass" (so called because it largely deter- 
mines the rate at which the system's frequency evolves, 
or "chirps"), / = O/tt = 2/ or b is the GW frequency, 
and $(t) = f 2nf(t')dt' = 2$ orb is the GW phase. We 
have chosen a coordinate system oriented such that the 
binary's orbit lies within the xj/-plane; this tensor will 
later be projected onto polarization basis tensors to con- 
struct the measured polarizations and h x - 

The rate at which the frequency changes due to the 
emission of gravitational radiation is given by (46| 



4f 

dt 



9(i 



11/3 



+ (47r-/3)(7rM/) 
/34103 13661 



59 



77 + — ?7 2 + a (ttM/) 4/3 
V 18144 2016 ' 18 ' J ' 

"(2.4) 

Notice that the chirp mass M. dominates the rate of 
change of /; the reduced mass \i and parameters (3 and 
a have an influence as well. The parameter j3 describes 
spin-orbit interactions and is given by 
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2 r 



/ 777, \ 2 




U3 (l?) 







L-S,- 



(2.5) 



The parameter a describes spin-spin interactions: 



48M(™ 2 m|) 



[721(L-Si)(L-S 2 )-247(Si-S 2 )] . (2.6) 
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Notice that j3 and a depend on the angles between the 
binary's angular momentum and the two spins. In pre- 
vious analyses which have neglected precession, (3 and a 
are constants; precession makes them time dependent. 
Using Eq. (|2.4[) , we can now integrate to find 



*(/) 



3 1 336 4 7 J ' 



+2 



> 

/3058673 5429 617 



V 1016064 1008 ' 144 



\4/3 



(2.7) 



The parameter t c formally defines the time at which / di- 
verges within the post-Newtonian framework. In reality, 
we expect finite size effects to significantly modify the bi- 
nary's evolution as the members come into contact. The 
system evolves so quickly as the bodies come together 
that t c is nonetheless a useful surrogate for a "time of co- 
alescence". Finally, the wave phase $(f) = J 2irf(t')dt' 
as a function of / is given by 



*(/) - *[*(/)] 
x 1 



1 

16 



(*Mf) 



-5/3 



5 /743 11 > , ,-, N2 /3 
3 + 



+5 



-(47r-/3)(7rM/) 
/ 3058673 5429 



V 1016064 1008 



'/ + ^ a) (ttM ff* 
(2. 



where $ c is the phase at time t c . The restricted PN 
waveform is then constructed by inserting (|2.8p into I 



B. Precession equations 

We next examine the effects of precession on the binary 
system. As discussed in the Introduction, spin-orbit and 
spin- spin interactions cause the black hole spins Si and 
S2 to precess. Precession occurs, at leading order 6 , on a 
time scale T plec cx r 5 / 2 at large separations [5(J. Since 
this is smaller than the inspiral time scale Ti nsp , we treat 
the total angular momentum J = L + S1 + S2 as constant 
over T prec . The orbital angular momentum L must then 



6 Several effects are built in to the precession equations discussed 
below, leading to precessions that occur on time scales scaling as 
r 5 / 2 and r 3 . Since we integrate these equations numerically, all 
of these effects are included in our analysis. For the purposes of 
this discussion, we subsume these effects into the leading-order 



precess to compensate for changes in Si and S2- Since 
T prcc is longer than the orbital time scale T or b , we use an 
orbit-averaged version of the precession equations [29j : 



S, = ^ 



2 + 

2 mi / 
is 2 -|(S 2 -L)L 

2 m 2 / 

-Si - -(Si-L)L 

2 2 V ' 



x Si 
x Si, 

x S 2 
x S 2 , 



(2.9) 



(2.10) 



where 7 r = M 1 / 3 /(7r/) 2 / 3 . These equations each have 
two pieces [12j. Consider the equation for Si. The first 
piece, which contains no S 2 dependence, is the spin-orbit 
term. This term, which comes in at 1PN order, is due 
to the geodetic precession of Si as hole 1 orbits in the 
spacetime generated by the mass of hole 2, and to the 
Lense-Thirring precession of Si in the gravitomagnetic 
field generated by the orbital motion of hole 2. The sec- 
ond piece is the spin-spin term, which enters at 1.5PN or- 
der. This term can be understood as the Lense-Thirring 
precession of Si in the gravitomagnetic field generated 
by the spin of hole 2. Note that the magnitudes of the 
spins do not change at this order; see [29| for more de- 
tails. From conservation of total angular momentum on 
short time scales, we have 



L = -(Si+S 2 ). 



(2.11) 



Over longer time scales, we must also consider the change 
in total angular momentum due to the radiation reaction, 
which is given by 



• 32 fi 2 (M 
5 r 



5/2 



(2.12) 



to lowest order. 

Considering only the spin-orbit terms and taking the 
limit S2 = or mi = rn 2 leads to a system whose preces- 
sion can described analytically; this is the "simple pre- 
cession" limit described in [2^| . Simple precession can be 
visualized as a rotation of L and S = Si + S 2 around 
the total angular momentum J. (Since inspiral shrinks 
J, the precession is actually around a slightly different 
direction Jo; see [29| for further discussion.) 

Since Vecchio restricts his analysis to mi = m 2 and 
does not include the spi n-sp in interaction, this limit is 
appropriate for his work [26]. As a consequence, Vecchio 
takes the quantities L • Sj, Si • S 2 , and /? to be constant. 



7 We use only the lowest-order Newtonian orbital separation in 
these equations. Including more terms would introduce higher- 
order effects into the precession. 
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(He does not include the spin-spin term a in the anal- 
ysis.) Here, we will study the impact of the full (albeit 
orbit-averaged) precession equations, including spin-spin 
terms, and include the impact of mass ratio. An ana- 
lytic description is not possible in this case, so we must 
integrate these equations numerically. The behavior is 
qualitatively similar to the simple precession case, but 
with significant quantitative differences. For example, 
P now oscillates around an average value. For unequal 
masses (say TO1/TO2 2), the difference due to precession 
can be substantial [47|. Such cases are also astrophysi- 
cally the most interesting — a mass ratio of roughly 10 is 
favored in binary black hole formation_scenarios arising 
from hierarchical structure formation (5lj |. 

The precession of the orbital plane causes a change in 
the orbital phase < & rb(i) [HI- Multiplying by a factor of 
2, the change in the wave phase is 

6 p <!>(t)=-J^ 5 p $(t')dt', (2.13) 

where 

2L n 

W) = - , f „. 2 (Lxn)-L, (2.14) 
1 — (L ■ n) z 

and n is the direction of the binary on the sky. At this 
point, we note that precession has several effects on the 
observed waveform. It changes the orbital phase, even 
at Newtonian order, by the amount (|2.13p . and it mod- 
ifies the functions (3 and a which appear in the post- 
Newtonian phase (|2.8p and time-frequency relation (|2.T|) . 
In the next section, we consider extrinsic effects on the 
waveform and find that precession of the orbital plane 
modifies them as well. 

C. Extrinsic effects 



These expressions were first discussed in the Introduction 
(albeit without the preccssional phase correction). Here 
Dl is the luminosity distance to the source. Notice that 
the weighting of the two polarizations depends upon the 
direction of the angular momentum vector relative to the 
sky position. 

We now consider the GW as measured by the detector. 
All of this analysis is done using the long wavelength 
(A 3> L, where L is the LISA arm length) approximation 
introduced by Cutler more details can be found there. 
This approximation is appropriate for our purposes since 
most of the signal accumulates at low frequencies where 
the wavelength is in fact greater than the arm length. 
The full LISA response function, including arm-length 
effects, is discussed in [52l l53j. 

LISA consists of three spacecraft arranged in an equi- 
lateral triangle, 5 x 10 6 km apart. The center of mass of 
the configuration orbits the sun 20° behind the Earth. 
The triangle is oriented at 60° to the ecliptic, so the or- 
bits of the individual spacecraft will all be in different 
planes. This causes the triangle to spin around itself as 
it orbits the sun. Following Cutler, we define a barred 
"barycenter" coordinate system (x, y, z), which is fixed in 
space with the xy-plane aligned with the ecliptic, and an 
unbarred "detector" coordinate system (x, y, z), which is 
attached to the detector. The z axis always points to- 
ward the sun, 60° away from vertical, while the x and y 
axes pinwheel around it. A particular binary will have 
fixed coordinates in the barycenter system, but its detec- 
tor coordinates will be time varying. 

The three arms act as a pair of two-arm detectors. We 
are first interested in the strain measured in detector I, 
that formed by arms 1 and 2: 



We have now constructed the intrinsic GWs emitted 
by a precessing binary in the restricted post-Newtonian 
approximation. The waveform measured by LISA will 
also include extrinsic effects due to the binary's location 
on the sky and the motion of the detector. 

We can write the wave as a combination of two or- 
thogonal polarizations propagating in the — n direction. 
Define p and q as axes orthogonal to n, with p = 
n x L/|n x L| and q = p x n. These are the principal 
axes for the wave; that is, they are defined so that the two 
polarizations are exactly 90° out of phase. The polariza- 
tion and 



tion basis tensors for these axes are = piPj 



P.'lj ■ 'LI'.- 



where 
h + (t) = 2 



h ij (t)=h + (t)H± + h x (t)H* j , (2.15) 
[1 + (L • n) 2 ] cos[$(i) + 6 p $(t)] , 



M 5/3 < f\2/3 , 



(2.16) 
(2.17) 



h(t) 



5L x (t) - 8L 2 {t) 



(2.18) 



where SLi(t) and 5L2(t) are the differences in length in 
arms 1 and 2 as the wave passes. L is the unperturbed 
length of the arms. Using the geometryof the detector 
and the equation of geodesic deviation [54]], we find 



hi(t) 



V3 
2 



(2.19) 



To obtain h xx and h yy for use in these equations, we 
must rotate the waveform from the principal axes into 
the detector frame. The result is that detector I measures 
both polarizations, modulated by the antenna pattern of 
that detector: 
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(On, $n, i/jn) cos[$(t) + (5 p $(t)] — 4(L-n)i ;i I >< (9n,4>n, V'jv) sin[$(t)+5 p $(t)]) . 

(2.20) 



Detector I acts like a "standard" 90° GW interferometer 
(e.g. LIGO), with the response scaled by \/3/2 (due to 
the 60° opening angle of the constellation) . The antenna 
pattern functions are given by 

F^(6 N ,(/) N ,ip N ) = i(l + cos 2 9 N ) cos20jv cos2ip N 



cos 9n sin 2<pN sin 2ip 



N 



(2.21) 



?f{dN,<t>N,i>N) = ^C 1 + COS"6' A r)cos2(/)ArSin2V'JV 



+ cos 6* n sin 207v cos 2?/>jv 



where On and </>j\r are the spherical angles for the binary's 
direction n in the (unbarred) detector frame and %pN is 
the polarization angle of the wave in that frame: 



tan ipN 



q • z L ■ z — (L • n) (z • n) 
P • z n • (L x z) 



(2.23) 



In order to use these expressions, we must relate the time- 
dependent angles in the unbarred detector frame, On, 4>n, 
and t/ijv, to quantities in the barred barycenter frame. 
Again, the details can be found in Cutler 

Similar expressions hold for detector II. Following Cut- 
ler, we construct the signal from detector II as 



hu(t) 



V3 



[h l (t)+2h lV (t)] , 



(2.24) 



where h\ is the signal from detector I (|2.18j) . and hw = 
(6Li2(t) — SLs(t))/L is the signal formed from the differ- 
ence in the lengths of arms 2 and 3. This choice makes 
the noise in detector I uncorrelated with the noise in de- 
tector II; we will exploit this property in Sec. HHI to treat 
detectors I and II as independent detectors. From <|2.24[) , 
we obtain 



hn(t) = ^ 



2 i^xy 



(2.25) 



The result is that detector II also behaves like a 90° in- 
terferometer (scaled by \/3/2), but rotated by 45° with 
respect to detector I. Thus the antenna patterns for de- 
tector II are 

F+(9 N , 4>n,4>n) = F+(6 N ,<j> N - tt/4, 4> n ) , (2.26) 
F* (On, <f>N, i> N ) = F{ (0 N , 4> n - tt/4, ■ (2.27) 

We now rewrite the waveform in terms of an amplitude 
and phase. Letting i = I, II label detector number, the 



waveform as measured by detector i is 



hi{t) = 2- ^U—A pol>i {t) C osm) 



(2.28) 



+ Ppol,i(*) + <fD(t) + S p <i>(t)} , 



where 



V3 r 



A poU (t) = ^-[(l + (L.n) 2 ) 2 ^ + (t) 2 +4(L.n) 2 F t X (tf]^ 



(2.22) j s "polarization amplitude" and 



¥> P ol,i(t) = tan" 



2(L-n)F i *(t) 
[l + (L-n)^+(t) 



(2.29) 



(2.30) 



is the "polarization phase" [f|. We have introduced the 
"Doppler phase" (fD{t), which arises from the detector's 
motion around the sun and is given by 



cp D (t) = 2wf(t)R (B smO N cos[$ D (t) 



(2.31) 



where $d(^) is the orbital phase of the detector and 
R® = 1 AU. 

Much of our analysis is done in the frequency domain. 
We define the Fourier transform of the signal as 



Hf) = 



h(t)dt . 



(2.32) 



To evaluate the Fourier transform, we make use of the 
stationary phase approximation [8j, l9j . This approxima- 
tion relies on the fact that the orbital time scale T or b 
is much shorter than the precession time scale T proc , as 
well as the inspiral time scale Ti nsp and detector orbital 
time scale Tn = 1 yr. The result thus differs from the 
true Fourier transform by terms of order T or b / Tp rec and 
Torb/Tlnsp 1^3 • ^ ne F° ur i er transform is thus likely to be 
inaccurate near the end of the inspiral, when all of these 
time scales become comparable. Using (|2.7p and (12.81) . 
we have 



= Vge— q^W*C/)]/ 7/6 (2.33) 



x e 



imf)-V vo x,iW))-<Pn[t{f))-8 p *[t{f)]) 
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where the phase ^(/) is given by 

*(/) = 2nft c $ c - - A + ±-^Mf)-^ 



4 128 
20 /743 11 



9 \336 

m/ 3058673 5429 

4 4tt - 3)(ttM f) + 10 1 77 

v MA I 1016064 1008 ' 



617 2 
+ 144" 



(2.34) 



In the work by Cutler 0, the separation of time scales 
that we used above leads to an interpretation of the polar- 
ization amplitude, polarization phase, and Doppler phase 
as modulations, in amplitude and phase, of an underly- 
ing carrier signal. These modulations make it possible to 
measure the sky position of the source, which also helps 
to measure the luminosity distance Dl H- With the 
addition of precession, the polarization amplitude and 
polarization phase include additional modulations which 
further improve the measurement of these parameters. 
In conjunction with the purely intrinsic effects of pre- 
cession, these effects also help us to better measure the 
masses and spins of the system. 



III. MEASUREMENT AND PARAMETER 
ESTIMATION WITH LISA 

A. Theory 

In the previous section, we constructed the expected 
form for the GW strain that LISA is being designed to 
measure. The signal Sj(t) as measured by detector i will 
of course also include noise rn(t): 



i {t) = h i {i)+n i {t). 



(3.1) 



The LISA noise spectrum is discussed in section MI Bl 
in this section, we discuss the theory of parameter es- 
timation with a noisy signal. First, consider only one 
detector. We assume that the noise is zero mean, wide- 
sense stationary, and Gaussian. Wide-sense stationary 
means that the autocovariance function 

K n {t,t') = {n{t)n{t')) - (n(t))(n(f)) (3.2) 

depends only on the time difference r = t — t'. (Through- 
out this section, quantities within angle brackets are en- 
semble averaged with respect to the noise distribution.) 
A process is Gaussian if every sample of the process can 
be described as a Gaussian random variable and all pos- 
sible sets of samples of the process are jointly Gaussian. 
However, the noise is colored, not white. A white noise 
process is defined to be a process which is uncorrelated 
with itself at different times; that is, its autocovariance is 
a delta function. Because the noise is colored, it has an 



interesting (nonflat) power spectral density (PSD), which 
is defined as the Fourier transform of the autocovariance 
function: 



Sn(f) 



dre 2 * lfT K n (r) . 



(3.3) 



The factor of 2 follows [8j ; we actually use the one-sided 
PSD. Since the noise is Gaussian, it is described entirely 
by its second moments. Therefore, we will only need the 
PSD, and not the full probability density function, to 
analyze the effect of the noise on the signal. 

Incidentally, it can be shown that wide sense station- 
arity implies that the Fourier transform of 77(7;) is a non- 
stationary white noise process in frequency: 



(n(f)n*(f')) = ls(f-f)S n (f)- 



(3.4) 



The Fourier components are thus independent Gaussian 
random variables. 

Now briefly consider both detectors. We explicitly con- 
structed the second detector (|2.24[) [with h(t) — > s(t)] so 
that the noise in it is uncorrelated with, and thus inde- 
pendent of, noise in the first detector. Thus we have 

(n l {f)n*{f)) = l -5 l AI-f')Sn{f). (3.5) 

The uncorrelated nature of these two noises will allow us 
to easily generalize discussion from one detector to the 
full two effective detector system. 

Let us write our generalized GW as h(6), where the 
components of the vector 9 represent the various param- 
eters on which the waveform depends. We now assume 
that a GW signal with particular parameters 9 is present 
in the data (i.e., "detection" has already occurred), and 
want to obtain estimates 9 of those source parameters. 
Finn [28[ shows that the probability for the noise to have 
some realization no(t) is given by 



p(n = 710) oc 



-(n |n )/2 



where the inner product used here is given by 



ra*(f)b(f) 



(a\b) = 4 Re / df 



Sn(f) 

a*(f)b(f)+a(f)b*(f) 



df- 



Sn(f) 



(3.6) 

(3.7) 
(3.8) 



This product is a natural one for the vector space of (fre- 
quency domain) signals a(f). (Note that this definition 
of the inner product differs from [28[ by a factor of 2.) 

Given a particular measured signal s(t), the probabil- 
ity that the GW parameters are given by 9 is the same 
as the probability that the noise takes the realization 
s-h(9): 



p{9\s) oc e 



-(h{0)-s\h(0)-s)/2 



(3.9) 
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where the constant of proportionality may include prior 
probability densities for the parameters 9. For simplicity 
we take these to be uniform. 

We can estimate the parameters 9 by the maximum 
likelihood (ML) method. This method involves finding 
the parameters 9 that maximize (|3.9[) . or alternatively, 
minimize (h(9) — s\h(9) — s), which can be considered a 
distance in signal space. A bank of template waveforms 
is correlated with the received signal and, assuming that 
any template produces a statistically significant correla- 
tion, the one with the highest correlation is the one with 
the ML parameters. The SNR for this signal is then given 

by a 



B. LISA detector and astrophysical noise 



(h(9)\h(9)Y/^ {h{9)\h{9))^ . 



(3.10) 



To quantify the errors in the ML estimate, we expand 
(|3.9p around the most likely values 9. We can then write 
the probability density as [1, Q : 



p(9\s) cx e- r <* sea5eb / 2 
where 58 a = 6 a - § a and 



dh 



dh 



(3.11) 



(3.12) 



evaluated at 9 = 9, is the Fisher information matrix. 
For small deviations from the ML estimate, the distribu- 
tion is Gaussian. This expression holds for large values 
of the SNR (|3.10p . It is worth emphasizing at this point 
that, in our evaluation of Eq. (|3.12| . most derivatives are 
taken numerically using finite differencing — the compli- 
cated nature of the signal (due to the inclusion of spin 
precession) makes it essentially impossible to evaluate all 
but a few of our derivatives analytically. This is another 
reason that the code we have developed for this analysis 
is substantially slower than those developed for analyses 
which do not include spin-precession physics. 

Now we return again to the two detector case. Using 
(|3.5p . we can write a total Fisher matrix as the sum of 
the individual Fisher matrices for each detector: 



rt0t T^I 
ab — 1 



ab 



ab • 



(3.13) 



The Fisher matrix is then inverted to produce the covari- 
ance matrix E ab = (r to J:) af) . The diagonal terms of the 
covariance matrix represent measurement errors: 



We turn now to a discussion of the noise we expect 
in LISA measurements. Our model for the instrumen- 
tal noise spectrum, S i £ st (f), is based on that described 
in Ref. [55j. (From now on, we use the notation Sh for 
strain noise instead of S n for general noise.) In particular, 
we use the online sensitivity curve generator provided by 
S. Larson, which implements the recipe of [5|| (see [56|). 
The output of Larson's webtool gives a sky averaged am- 
plitude sensitivity curve, /iLarson- To convert to the noise 
we need for our analysis, we square this amplitude and 
insert two numerical factors: 



canst 



m- 1 r^h 

(J) — g x I ""^""•Larson 




(3.16) 



The factor of 1/5 accounts for the averaging of the an- 
tenna pattern functions over all sky positions and source 
orientations. This factor is only correct for measuring 
radiation with wavelength A 3> L (where L is the LISA 
arm length). As a consequence, our instrumental noise 
will be inaccurate at high frequencies. This will have lit- 
tle impact on our analysis since, as already argued, the 
signal from merging binary black holes accumulates at 
low frequencies. 

The factor V3/2 arises due to the 60° opening angle 
of the interferometer arms; we have already accounted 
for this factor in our discussion of the interferometer's 
interaction with a GW [cf. Eqs. (|2~T§)) and (|235]) ]. The 
numerical factor 3/20 has been the source of some confu- 
sion; Berti, Buonanno and Will very nicely straightened 
this out. See Sec. IIC of Ref. [7j for further discussion of 
these factors. 

Besides purely instrumental noise, LISA data will 
contain "noise" from a background of confused binary 
sources 8 , mostly white dwarf binaries. An isotropic back- 
ground of indistinguishable sources can be represented as 
noise with spectral density [27l | 



neon 



(/) = — r-Vwco, 

07T 



(3.17) 



where p c = 3Hq/8it is the critical energy density to close 
the universe and flaw = {f / Pcjdpcw / df is the energy 
density in GWs relative to p c per logarithmic frequency 
interval. Using this form and the results of Farmer and 
Phinney (57l | , we model the confusion noise due to extra- 
galactic binary sources by 



Ma - V/PW = VE= ■ (3.14) ^ n _ Ao ^ n _ 47 " V» ^ , 

The off-diagonal terms can be expressed as correlation 
coefficients, ranging from —1 to 1: 



Sr sal (/) = 4.2 x 10- 47 -4- Hz" 1 . (3.18) 

1Hz/ 



c = \a A ah = ~1^=^= ■ (3.15) 8 While surely noise when studying cosmological black holes, this 

Ap AW vS aa S^ background is signal to those interested in stellar populations. 
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From Nelemans et al. [58| . we take the galactic white 
dwarf confusion noise to be 

Sf(f) = 2.1 x 1(T 45 U^j V3 Hz" 1 . (3.19) 

The combined instrumental and galactic confusion noise 
is given by [13] 

5 in st+g al (/) = min[ ^ns t(/)/exp( _ KT -l siond7V/d/); 

sr\f)+sf(f)}. 

(3.20) 

The choice taken in Eq. (|3.20|) reflects the fact that, at 
sufficiently high frequency, the number of binaries per 
bin should be small enough that they are no longer truly 
confused and can be subtracted from the data stream 
(at least partially). The factor KT~^ ssion dN / df) is 
the fraction of "uncorrupted" frequency bins. We choose 
k = 4.5 [59l ] . Tmission is the mission duration (which we 
take to be three years), and 

^= 2)<1 0-(i^)" /3 H z - (3.21) 

is the number density of galactic binaries per unit fre- 
quency ((J. 

Finally, the total noise is given by 

S h (f) = Sr t+ ^(f) + S c h ^ l (f). (3.22) 

IV. RESULTS 
A. Procedure 

1. Parameter space 

Seventeen parameters describe the most general binary 
black hole inspiral waveform [26|. Two of these are the 
orbital eccentricity and the orientation of the orbital el- 
lipse; since we only consider circular orbits, we can ignore 
these two. The other 15 parameters are all necessary to 
describe the full post-Newtonian waveform with preces- 
sion effects that we described in section [TT1 

We divide this set into intrinsic and extrinsic param- 
eters. In our labeling system, intrinsic parameters are 
those which label properties intrinsic to the binary it- 
self; extrinsic parameters label properties which depend 
upon the position and placement of the binary relative 
to the observer. One can regard intrinsic parameters as 
describing the physics or astrophysics of the binary sys- 
tem, and extrinsic parameters as describing the binary's 
astronomical properties. 

The intrinsic parameters we use are In mi; In 7712; 
/2l(0) = cos[#l(0)] and </>l(0), the initial direction of the 
orbital angular momentum; p,Si(0) = cosf^ (0)], 0s 1 (O), 
/2s 2 (0) = cos[(?s 2 (0)], and <^s 2 (0), the initial directions 



of the spins; xi and X21 the dimensionless spin param- 
eters; t c , the time at coalescence; and $ c , the phase at 
coalescence. (Note that t c and $ c could very well be con- 
sidered extrinsic in our labeling system, since they just 
label the system's state at some particular time. At any 
rate, neither t c nor <i> c is of much physical interest, so 
their categorization isn't too important.) Our extrinsic 
parameters are £Ln — cos6*jv and 4>n, the sky position 
in barycenter coordinates; and lnD^, the luminosity dis- 
tance to the binary. All of these parameters must be 
fit in a measurement and thus must be included in our 
Fisher matrix analysis. We are not necessarily interested 
in all of them, however. In particular, we will focus on 
the masses, the dimensionless spin parameters, the sky 
position, and the luminosity distance. 

It is worth noting that this choice of parameters is not 
the same as that used in analyses which neglect preces- 
sion. In that case, the direction of the angular momen- 
tum L is constant, and so the system's orientation is con- 
stant and fully described using two angles (e.g., fir, and 
4>l)- Including precession, L is no longer constant, but 
evolves according to (|2.11[) . The solution to this differen- 
tial equation requires two initial conditions, for instance, 
/2l(0) and 0^(0), which can be used as parameters of the 
system. Since these initial conditions are taken at the 
(somewhat arbitrary) starting point of our calculations, 
they do not hold much physical interest (though they 
must be fit for and thus included in our Fisher matrix). 

Previous analyses, including the precursor to this work 
@ , have used (5 (|2.5|) and a (|2.6|) as parameters — these 
are constants when precession is neglected. They are also 
the only combinations of the spin magnitudes and spin 
angles that enter into the expression for the waveform. 
Boiling the six numbers which characterize Si and S2 
down to two greatly simplifies the parameter space, but 
also restricts us from being able to measure, for exam- 
ple, the black holes' spin magnitudes. When precession 
is included, (3 and a are no longer constants. In addi- 
tion, they no longer fully characterize the signal, since 
the precession equations (|2.9[) . (|2.10p . and (|2.11[) depend 
on all of the components of the spins. We thus need six 
spin-related parameters to fully describe the signal: the 
magnitudes of the spins and their orientations at some 
initial time. The orientations are again uninteresting, 
but the fact that we can measure the magnitudes of the 
spins and quantify their errors is quite interesting and 
new to this analysis. 

Finally, we break from tradition and use In mi and 
lnm,2 to parameterize our masses rather than ln.M and 
ln^. The chirp mass and reduced mass have been used 
in most previous work because of their appearance in the 
waveform phase ^(Z). However, the precession equa- 
tions, as well as the spin parameters (3 and er, depend 
on the individual masses of the black holes. It is a sim- 
ple matter in principle to just solve for mi,2(.M,/i) and 
substitute into the precession equations. Unfortunately, 
the Jacobian of the transformation between (M , /i) and 
(mi, m.2) is singular when mi = m,2, leading to problems 



13 



in evaluating the Fisher matrix. 

These problems can be illustrated analytically. Con- 
sider how derivatives of some function /(mi, m^) with 
respect to M. behave: 



df _ \ - df(m 1 ,m 2 ) dmj{M,n) 
dM ~ dm, dM ' 1 ' ' 

i—l 

When mi — m 2 , the final derivative diverges — a be- 
havior that we have seen numerically. The Fisher in- 
formation is infinite, and the Gaussian approximation 
breaks down; the same problem occurs for fi. Thus, we 
argue that, when precession is included, M and /j, are no 
longer a good choice of parameters to describe the sys- 
tem. Since we are still interested in the errors in hiM 
and In fj, (which are determined to higher accuracy than 
the individual masses), we convert using the propagation 
of errors formulas 




For unequal masses, we find that computing errors in 
mi and 7712 and then converting gives the same result as 
simply computing errors in M and /1 directly. We do 
not find good agreement in the equal mass case; for the 
reasons discussed above, however, we do not trust the 
(M,n) parameterization in this case. At any rate, the 
case mi = m,2 is quite implausible in nature, so this is 
almost certainly a moot point as far as real measurements 
are concerned 9 . We note that Vecchio [26[, for simplicity, 
considers the equal mass case exclusively but does not 
report any anomalous behavior such as we have seen. 
We are puzzled about this discrepancy. 



It is worth noting that even a slight mass difference (a few per- 
cent) is sufficient for the two approaches to match. 



2. Calculations 

The code we use to calculate parameter measurement 
errors is based on that used in @. It is written in C++ 
using several routines taken, sometimes with slight mod- 
ification, from As in jfjj], we perform Monte Carlo 
simulations in which we specify rest-frame masses and 
redshift and then randomly choose sky position, initial 
angular momentum and spin directions, spin magnitudes, 
and time of coalescence within the three- year LISA mis- 
sion window. We specify spin magnitudes for some stud- 
ies as well. 

The primary function of the code is the calculation of 
the full gravitational waveform, including precessional ef- 
fects. In order to effectively use the formulas of section 
IIIIl we take the wave frequency / as the independent 
variable. The elapsed time is related to the frequency 
using (|2.7p . The calculation is started when the wave- 
form enters LISA's band (taken to be f m - ln = 3 x 1CP 5 
Hz throughout this paper) or when the LISA mission 
begins, whichever is later. (By treating the time of coa- 
lescence as a Monte Carlo variable, some signals will be 
partially cut off because they are already in band when 
LISA begins observations.) The calculation proceeds un- 
til the binary reaches the Schwarzschild innermost sta- 
ble circular orbit (ISCO) at orbital separation r = 6M. 
Though perhaps a somewhat crude choice, we use this 
criterion for simplicity. Choosing slightly different cut- 
off radii does not change our results very much; at any 
rate, the post-Newtonian phase formula and precession 
equations we use in this regime are unlikely to be very 
accurate. The frequency at this point, which we call the 
"merge frequency," can be found using (|2.1|) for r = 6M 
(plus / = n/ir). For simplicity, we just use the lowest 
order term: / = (Af/r 3 ) 1 / 2 /n. Any error due to this 
approximation is no doubt unimportant compared to the 
arbitrary selection of r = 6M for the transition. 

Once the frequency range has been determined, the 
true work begins. We integrate the precession equations 
([23]) . (|2~TU1) . and (j2~TTj) using a Runge-Kutta routine to 
find the values of L, Si, and S 2 over the duration of the 
signal. The routine is a fifth-order adaptive step algo- 
rithm in the frequency domain. At each frequency, the 
code takes the results for the three orbital angular mo- 
mentum components and six spin components and uses 
them to calculate fl^, 4>l, (3, and a. It also computes the 
integrated correction to the phase using the derivative 

(EH. 

As already discussed, our derivatives are taken numer- 
ically rather than analytically. We therefore must do 
the integration described above a total of 21 times: once 
for the given values of the parameters, and twice more 
for small shifts in each parameter which requires a nu- 
merical derivative. This repetition slows the code quite 
drastically compared to its earlier incarnation — an un- 
fortunate but unavoidable cost. 

Once all of the necessary integrations are complete, 
the SNR (j3"+U|) and the Fisher matrix (j3"T2"|) can be cal- 
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culated for each of the two effective detectors of LISA 
using the noise Sh(f) (see Sec. IIIIB[) . Some previous 
work 5, 7] investigated parameter estimation using the 
signal from only one synthesized detector; we will always 
assume that both are operational. It would be interesting 
to see how measurement degradation due to only having 
a single operating detector can be ameliorated by includ- 
ing precession effects. 

At this stage, the necessary integrals are performed us- 
ing Curtis-Clenshaw quadrature, which depends on the 
decomposition of the integrand into Chebyshev polyno- 
mials [60|| . This method keeps the code reasonably fast 
even with the addition of the Runge-Kutta routine. At 
each step of the integration, the integrator uses the val- 
ues that were calculated using that Runge-Kutta routine 
to evaluate the waveform and/or its appropriate deriva- 
tives. The derivatives are calculated using 

d9~ A9 ' [ ' ' 

For all parameters, we use A9 = 10~ 5 9. We invert the 
Fisher matrix using LU decomposition to produce the 
covariance matrix [60(. In "poor" cases, (e.g., high mass 
binaries at large redshift) the Fisher matrix can be nearly 
singular, with a large condition number 10 . In such a 
case, the covariance matrix produced by the code may 
not be the true inverse of the Fisher matrix (and may 
not even be positive definite). This problem is largely 
ameliorated by representing our numerical data in long 
double format — this improves (relative to type double) 
matrix inverses in many "bad" cases but leaves all other 
cases essentially unchanged. 

It is worth noting that the bad cases are typically ones 
in which the binary executes very few orbits over the 
course of the measurement. We are confident in our re- 
sults for all cases in which the number of measured orbits, 
iV or b = $ orb/27r, is greater than ~ 10 — 20. When the 
number of orbits is small (and the condition number is 
concomitantly high), the errors are so large that they 
are basically meaningless. In such a case, measurement 
would not determine the system's characteristics in any 
meaningful sense. 



B. Black hole masses and spins 

Representative examples of our results are shown in 
Figures [3 and [3] These histograms show the spread of er- 
rors in Ai and fj, for a sample of 10 4 binaries at z = 1 with 
rest frame masses mi = 10 6 M and m-i = 3 x 10 5 M©. 



The "condition number" is the ratio of the largest eigenvalue of a 
matrix to the smallest. A rule of thumb is that matrix inversion 
breaks down when the logarithm of the condition number of a 
matrix exceeds the number of digits of accuracy in the matrix 
elements (see, e.g., discussion in [6011 ). 



Each figure compares the results of the new code to 
those of the original code of ||, which neglects preces- 
sion. (This code has been updated to reflect up-to-date 
models for LISA noise; some minor coding errors have 
also been corrected.) Clearly, including spin precession 
leads to a significant improvement in the measurement of 
these mass parameters. The reduced mass (j,, in partic- 
ular, is improved. This is because the time variation of 
j3 and a breaks a near degeneracy between those terms 
and /j, in the post-Newtonian phase (|2.34[t . The masses 
also control the precession rate, as seen in (|2.9[) . (|2.10[) . 
and (|2.1ip . (Recall that, in those equations, Si = Xi m j-) 
This means that they now influence the polarization am- 
plitude and polarization phase; they do not influence 
those quantities when precession is neglected. These 
precession-induced influences on the waveform make it 
possible to determine the masses even more accurately 
than before. 
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FIG. 2: Distribution of errors in chirp mass M for 10 4 bi- 
naries with mi = 1O 6 M and mi = 3 x 10 5 M Q at z = 1. 
The dashed line is the precession-free calculation; the solid 
line includes precession. Precession reduces the measurement 
error by about an order of magnitude. 

As discussed earlier, we have found the masses mi and 
TO2 to be more useful parameters than M. and [i when 
precession is included. Figure 2] shows the error in mea- 
surements of the individual masses for our example sys- 
tem. While these masses are measured quite accurately, 
they are not measured as accurately as M. and \i. This 
reflects the fact that, even though the individual masses 
play a role in the precession, the other parts of the wave- 
form depend explicitly on the combinations M. and /z. 
Notice also that the smaller mass is typically determined 
a bit better than the larger one, though the difference is 
not large. 

Precession makes it possible to determine the spins of 
the binary's members. Figure [5] shows the error in mea- 
surements of the two dimensionless spin parameters xi 
and X2- We see that x is generally determined very well: 
Taking a typical spin parameter to be about 0.5 (recall 
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FIG. 3: Distribution of errors in reduced mass for 10 4 bi- 
naries with mi = 1O 6 M and m,2 = 3 x 10 5 M Q at z = 1. 
The dashed line is the precession-free calculation; the solid 
line includes precession. Precession has an enormous effect 
on the reduced mass, which was previously highly correlated 
with the parameters (3 and a. 
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FIG. 5: Distribution of errors in dimensionless spin parame- 
ters xi (solid line) and \2 (dashed line) for 10 4 binaries with 
mi = 10 6 Mq and m2 = 3 x 10 5 A/ Q at z = 1. In each binary, 
the spin values are randomly selected between and 1. The 
higher mass then has, on average, higher total spin and more 
effect on the precession. 
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Anii/mi 

FIG. 4: Distribution of errors in individual hole masses for 
10 4 binaries at z = 1. The solid line is mi = 1O 6 M , while 
the dashed line is rriz = 3 x 10 5 A'/ Q . The individual masses 
are not determined as well as M and ft, but they are better 
behaved parameters when precession is introduced. 



we randomly choose \ between and 1), the bulk of this 
distribution corresponds to errors of a bit less than a 
percent. For this entirely random distribution of the 
dimensionless spin parameter of the larger hole tends to 
be better determined than that of the smaller hole. This 
appears to be a simple consequence of the fact that black 
hole spin scales as mass squared (Si = Xi m i)j an d larger 
spin has more of an impact on the waveform. 

Next, we examine how well spin is measured as a func- 
tion of spin magnitude. Figure [5] shows the error in 
Xi for the same system as in Fig. [5J except that we 



set Xi = X2 = 0.9 (solid line) and Xi — X2 = 0.1 
(dashed line), rather than randomly distributing their 
values. This allows us to more accurately assess how well 
spin is determined as a function of its value, as well as to 
more accurately determine the percent error we expect 
in these measurements. For Xi — X2 — 0.1, the error is 
almost 10%, while for xi — X2 = 0.9, the error is closer 
to 0.1%. This is a considerable difference and is easily 
ascribed to the fact that rapid spin has a much stronger 
impact on the waveform. 

Table|T]shows the median errors in intrinsic parameters 
for different mass ratios at z = 1. We continue to include 
the errors in M. and ijl for comparison with the precession- 
free case, but only in binaries of unequal mass where the 
Gaussian approximation is well defined. Examining the 
table, we see some interesting features. The errors, in 
general, are worse for higher mass binaries, which spend 
less time in the LISA band. At mi = m% — 10 7 M Q , the 
mass errors jump to nearly 10%, compared to tenths of a 
percent at the next lower mass combination. In addition, 
the spin determination becomes very unreliable. Mass 
ratio also has an important effect on the results. Tak- 
ing into account the general trend caused by total mass, 
we see that unequal mass ratios generally produce bet- 
ter results. This is good news for eventual measurements 
of astrophysical systems, since merger tree calculations 
show that binaries are most likely to have mass ratios 
of about 10 [5l|. To understand the mass ratio depen- 
dence, we again turn to the precession equations (|2.9[) . 
(|2.10[) . and (|2.11[) . For unequal mass ratios, the geodetic 
spin-orbit and spin-spin terms will cause the two spins 
to precess at different rates, creating richer features in 
the signal than for equal mass ratios. This illustrates the 
importance of effects beyond the "simple precession" of 
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FIG. 6: Distribution of errors in dimensionless spin parameter 
Xi for 10 4 binaries with mi = 10 6 A/ Q and mi = 3 X 10 5 M Q 
at z = 1. Here, spin magnitudes have been set to a specified 
value — low spin, Xi — X2 — 0-1 (dashed line), and high spin, 
Xi — X2 = 0.9 (solid line). Since greater spin more strongly 
impacts the waveform, the high spin case is measured more 
accurately. 



Hallal- We also see that the trends of Figs. H and M hold 
for each unequal mass binary in the table. That is, the 
mass of the smaller hole is determined better than the 
mass of the larger hole, but the spin of the larger hole is 
determined better than the spin of the smaller hole. 

Tables [TT] and IIIII show the same results for z = 3 and 
z — 5, respectively. The trends we see at z = 1 largely 
continue at these redshifts. In general, the errors get 



worse at higher redshift as the signal amplitude degrades 
and more of the signal is rcdshiftcd out of band. It is 
worth noting that the change is generally greater from 
z = 1 to z = 3 than from z = 3 to z = 5. This 
effect was also seen by Berti, Buonanno, and Will Q 
and can be explained by considering the redshift depen- 
dence of the wave amplitude. Neglecting all the angular 
factors and remembering to redshift quantities with the 
dimensions of time, we find that the amplitude scales 
like (1 + z)/D L (z) = l/D M (z), where D M {z) is the 
proper motion distance. This distance measure varies 
more strongly with z at low redshift than at high red- 
shift. (See [61[ for a plot of Dm(z).) Consequently, when 
moving from z = 1 to z = 3, the amplitude, and thus the 
SNR, decreases more than when moving from z — 3 to 
z = 5. For lower mass binaries, this amplitude decrease 
plays a bigger role in the loss of SNR than does redshift- 
ing the spectrum to lower frequency; most of the SNR is 
accumulated late in the inspiral, where the orbits are in 
a relatively flat region of the sensitivity curve. 

By contrast, for the highest mass binaries, redshifting 
of the spectrum can have a dramatic effect. So much of 
their signal is moved out of band that LISA may measure 
their waves for only a very short time. As such, measure- 
ment may not provide sufficient information to constrain 
15 parameters. This is reflected in the high condition 
numbers associated with such cases. Their Fisher matri- 
ces are thus nearly singular, and their inverses are un- 
trustworthy. In fact, measurement error in these binaries 
actually degrades when precession is included. The time 
in band is too short for precession effects to accumulate. 
They do not aid parameter estimation; instead, the need 
to fit extra parameters causes errors to be worse. 



C. Sky position and distance to source 

We now focus on extrinsic parameters, the sky posi- 
tion and the luminosity distance to the source. We find 
that the determination of these parameters is likewise 
improved when precession physics is taken into account, 
though not as strongly as for intrinsic parameters. This 
might be expected, since precession is an intrinsic effect 
local to the binary and has no direct dependence on these 
extrinsic parameters. Precession's impact on the extrin- 
sic parameters is somewhat more indirect — it largely 
improves their determination by reducing the (otherwise 
quite strong) correlation between sky position and the or- 
bital angular momentum direction L and between these 
angles and the source's luminosity distance. 

In our analysis, a binary's position on the sky is char- 
acterized by the two parameters = (/2jy = cosd^, 4>n). 
Consider the subspace containing just these two param- 
eters. The Fisher matrix is a 2 x 2 matrix which char- 



acterizes the probability density for the true parameters 
given a measured signal [see p. lip ]: 

p(9\s) cx exp (~r ab 50 a 6^'\ . (4.5) 

To accurately describe the error ellipse on the sky, we 
need to manipulate the right hand side of (|4.5|) in sev- 
eral ways. First, we change coordinates from p,^f to 
On, while leaving (f>N alone; schematically, we can rep- 
resent this as 9 a — > a . This transformation gives 
5@n = (d@N /dp.N)SpN — —^P-n/ sinfljy 

Next, we need to redefine what we mean by "error" 
in order to make the results more relevant to observa- 
tions. To do so, we define the "proper" angular errors 
S0 N = 56n and 84> P N — sitiOnSiPn ■ The proper angu- 
lar errors are just the normal coordinate errors rescaled 
by the metric of the sphere to correctly account for the 
proper size of a segment 5(f) n at different On- Substitut- 
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mi (M ) 


m 2 (M Q ) 


Ami/mi 


Am2/m2 




AX2 


AM/M 
(no precession) 


AM/M 
(precession) 


Afi/fi 
(no precession) 


Afi/fj, 
(precession) 


10 5 


10 5 


0.000783 


0.000782 


0.00415 


0.00414 










3 x 10 5 


10 5 


0.000667 


0.000541 


0.00157 


0.00306 


5.92 x 10~ 5 


5.51 x 10~ 6 


0.0114 


0.000239 


3 x 10 5 


3 x 10 5 


0.00109 


0.00109 


0.00539 


0.00536 










10 6 


10 5 


0.000629 


0.000440 


0.00102 


0.00440 


0.000156 


1.18 x 10 -5 


0.0180 


0.000343 


10 6 


3 x 10 5 


0.00111 


0.000882 


0.00256 


0.00499 


0.000170 


1.19 x 10" 5 


0.0274 


0.000423 


10 6 


10 6 


0.00195 


0.00195 


0.00902 


0.00897 










3 x 10 6 


3 x 10 5 


0.000988 


0.000691 


0.00137 


0.00563 


0.000583 


2.53 x 10" 5 


0.0550 


0.000539 


3 x 10 6 


10 6 


0.00238 


0.00192 


0.00380 


0.00674 


0.00117 


4.19 x 10~ 5 


0.135 


0.000849 


3 x 10 6 


3 x 10 6 


0.00584 


0.00582 


0.0271 


0.0275 










10 7 


10 6 


0.00239 


0.00177 


0.00233 


0.0122 


0.00770 


0.000174 


0.469 


0.00140 


10 7 


3 x 10 6 


0.00814 


0.00671 


0.00829 


0.0159 


0.00851 


0.000436 


0.607 


0.00332 


10 7 


10 7 


0.0804 


0.0802 


0.492 


0.493 











TABLE I: Median errors in intrinsic quantities for 10 4 binaries of various masses at z = 1, including comparisons with the "no 
precession" case where possible. We have omitted the errors in chirp mass and reduced mass for equal mass binaries because 
that parameterization of the waveform fails the Gaussian approximation at those points. 



mi (M s ) 


m 2 (M Q ) 


Ami/mi 


Am 2 /m 2 


Axi 


AX2 


AM/M 
(no precession) 


AM/M 
(precession) 


An/n 
(no precession) 


Afi/n 
(precession) 


10 5 


10 5 


0.00362 


0.00362 


0.0187 


0.0185 










3 x 10 5 


10 5 


0.00363 


0.00294 


0.00879 


0.0171 


0.000406 


3.31 x 10~ 5 


0.0715 


0.00130 


3 x 10 5 


3 x 10 s 


0.00569 


0.00569 


0.0271 


0.0269 










10 6 


10 5 


0.00330 


0.00231 


0.00498 


0.0208 


0.00120 


7.09 x 10~ 5 


0.128 


0.00180 


10 6 


3 x 10 s 


0.00648 


0.00517 


0.0120 


0.0229 


0.00174 


9.17 x 10~ 5 


0.228 


0.00248 


10 6 


10 6 


0.0138 


0.0139 


0.0627 


0.0630 










3 x 10 6 


3 x 10 5 


0.00569 


0.00402 


0.00664 


0.0287 


0.00633 


0.000241 


0.456 


0.00314 


3 x 10 6 


10 6 


0.0181 


0.0148 


0.0223 


0.0386 


0.00708 


0.000554 


0.596 


0.00658 


3 x 10 6 


3 x 10 6 


0.0744 


0.0737 


0.412 


0.415 










10 7 


10 6 


0.0301 


0.0283 


0.0256 


0.177 


0.0189 


0.00506 


0.690 


0.0231 


10 7 


3 x 10 6 


0.434 


0.359 


0.282 


0.448 


0.0182 


0.0428 


0.643 


0.180 


10 7 


10 7 


12.1 


12.0 


62.2 


61.5 











TABLE II: Median errors in intrinsic quantities for 10 4 binaries of various masses at z = 3. 



ing all of these changes into (|4.5p . we obtain 



K0|s)(xexp(--r^,5« 



(4.6) 



where we have defined a proper Fisher matrix for the 
parameters (On^n)- In terms of the original Fisher ma- 
trix, the elements are 



T P 

On 


= = sin 2 


^N^fJ, N fJ.N I 


(4.7) 


T P 

0n<i 


- = r p - 


— r 

?JV ^N<t>N > 


(4.8) 


T l 

<PN<i 


- = CSC 2 

>N 




(4.9) 



Finally, we diagonalize the Fisher matrix by rotating our 
parameterization, 9 a — > 6 a , such that the probability 
(14. 6p becomes 



p(6\s) tx exp 



2\2 



(SOD 2 , (591) 



2«) 2 2(af) 2 



In these coordinates, the covariance matrix is 



K) 2 







P\2 



(4.10) 



(4.11) 



Following Cutler |5| , we define the error ellipse such that 
the probability that the source lies outside the error el- 
lipse is e . The semiaxes of the error ellipse are given 
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mi (M Q ) 


m 2 (M Q ) 


Ami /mi 


Ami /mi 




AX2 


AM/M 
(no precession) 


AM/M 
(precession) 


Afi/fj, 
(no precession) 


A/i//x 
(precession) 


10 5 


10 s 


0.00791 


0.00792 


0.0392 


0.0389 










3 x 10 5 


10 5 


0.00811 


0.00658 


0.0193 


0.0359 


0.00103 


8.00 x 10" 5 


0.172 


0.00290 


3 x 10 5 


3 x 10 5 


0.0134 


0.0134 


0.0615 


0.0616 










10 6 


10 5 


0.00718 


0.00502 


0.00993 


0.0409 


0.00326 


0.000184 


0.305 


0.00391 


10 6 


3 x 10 5 


0.0156 


0.0124 


0.0249 


0.0460 


0.00427 


0.000289 


0.469 


0.00596 


10 6 


10 e 


0.0424 


0.0423 


0.197 


0.200 










3 x 10 6 


3 x 10 5 


0.0161 


0.0117 


0.0158 


0.0808 


0.0115 


0.00103 


0.643 


0.00922 


3 x 10 6 


10 6 


0.0576 


0.0475 


0.0606 


0.107 


0.0108 


0.00265 


0.635 


0.0214 


3 x 10 6 


3 x 10 6 


0.396 


0.391 


2.43 


2.44 










10 7 


10 e 


0.279 


0.282 


0.208 


1.41 


0.0374 


0.0640 


0.704 


0.232 


10 7 


3 x 10 6 


10.1 


8.41 


6.10 


7.61 


0.106 


1.11 


0.769 


4.28 


10 7 


10 7 


2280 


2290 


10300 


9900 











TABLE III: Median errors in intrinsic quantities for 10 4 binaries of various masses at z = 5. The results for the highest masses 
are meaningless — the parameters are completely undetermined. 



by ^/2(cr? g) 2 . These quantities follow from the eigenval- 
ues of the covariance matrix (|4.11[) ; since eigenvalues are 
invariant under rotation, we can calculate them before 
performing the rotation. In terms of our original covari- 
ance matrix, the major axis 2a and minor axis 2b of the 
ellipse are given by 



esc 2 fl/vS^™ + sin 2 6 N E* N * N ± 



(esc 2 6>jv£>w - sin 2 6*ArE^«*«) 2 + 4(I>« < ?«) : 



1/2 



(4.12) 

taking the plus and minus for major and minor axes, 
respectively. We also find the area of the error ellipse: 



Att N = nab = 2?ry - (^n<Pn)2 _ (4.13) 

Many previous analyses have reported the ellipse's area 
AD, n or y/Ailjy, the side of a square of equivalent area, 
as the sky position error |§, 0, E3|- Information about 
the ellipse's shape, crucial input to coordinating GW ob- 
servations with telescopes, is not included in such a mea- 
sure. By examining both 2a and 2b, this information is 
restored. Figure [7] shows the major axis of the error el- 
lipse 2a for both the original code, with no precession, 
and the code including precession effects. Figure [5] shows 
the same for the minor axis 2b. (Note that these figures 
cannot tell us which major axis is associated with which 
minor axis; that information is lost in the construction 
of the histograms.) 

Compared to the code which does include precession 
physics, the median of both distributions is reduced by 
about half an order of magnitude. The minor axis dis- 
tribution also shows a long tail of very small errors. In 
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FIG. 7: Distribution of the major axis of the sky position 
error ellipse, 2a, for 10 4 binaries with mi = 1O 6 M0 and 7712 = 
3 x 1O 5 M0 at z — 1. The dashed line is the precession- free 
calculation; the solid line includes precession. Sky position, 
as an extrinsic parameter, is improved somewhat indirectly 
by precession; therefore, the improvement is less than for the 



those cases, the position would be very well-constrained 
in one direction. 

Finally, we examine how well distance to the binary is 
determined. Figure [9] compares ADl/Dl both with and 
without precession physics taken into account. For this 
case, the distance error improves by about a factor of 3. 

Table HVl shows the median extrinsic errors for binaries 
of different mass. For comparison purposes, we include 
results that neglect spin precession. Binaries with the 
best determined parameters at this redshift have total 
mass several x 10 5 M fl < M tot < several x 10 6 M n - 
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FIG. 8: Distribution of the minor axis of the sky position 
error ellipse, 26, for 10 4 binaries with mi = 1O 6 M0 and m.2 = 
3 x 1O 5 M0 at z = 1. The dashed line is the precession- free 
calculation; the solid line includes precession. 

smaller binaries are not quite determined so well due to 
the weakness of their signal, while larger ones are not 
determined so well because they radiate fewer cycles in 
band. We also see again that unequal mass binaries give 
better results than equal mass binaries due to the im- 
pact of mass ratio on precession effects. Overall, we find 
that the major axis of the error ellipse is on the order of 
a few x 10 arcminutes, while the minor axis is a factor of 
2 — 4 smaller. This represents an improvement over the 
"no precession" case by a factor ~ 2 — 7 for the major axis 
and a factor ~ 2 — 10 for the minor axis. The distance 
errors are on the order of 0.2% — 0.7% for most masses, 
a factor of ~ 2 — 7 improvement. 



Tables [V] and IVII show the same results for higher red- 
shift. We see the same trends as at z = 1, but with some 
degradation in numerical value. The sky position errors 
reach a few degrees in the major axis and several tens of 
arcminutes up to a degree or two in the minor axis. The 
distance errors are on the order of 1 to several percent 
for most masses. At the highest masses, we again see 
that these parameters are essentially undetermined and 
that precession makes things worse by requiring extra 
parameters to be fit. 



1200 




FIG. 9: Distribution of errors in the luminosity distance for 
10 4 binaries with mi = 1O 6 M and m 2 = 3 x 10 5 M o at 
z — 1. The dashed line is the precession- free calculation; the 
solid line includes precession. 



V. SUMMARY AND CONCLUSIONS 

The general relativistic precession of black holes in bi- 
nary systems can have a strong influence on the binary's 
dynamics @, [5(| and thus upon the GWs that it gen- 
erates. It has been known for some time that it will 
be necessary to take these dynamics into account in or- 
der to detect these black holes in noisy detector data 
EUSmSlH,!!!!!!; clearly, taking these dynamics 
into account will be just as (if not more) important for 
the complementary problem of determining the parame- 
ters which characterize a detected system. Vecchio [26[ 
first demonstrated that, by taking into account preces- 
sion physics, quite a few near degeneracies among binary 
source parameters can be broken, making our estimates 
for how accurately they can be determined more opti- 
mistic. This analysis largely confirms and extends Vec- 



chio's pioneering work. By taking the equations of mo- 
tion to higher order to include spin-spin couplings, and 
by surveying measurement accuracy as a function of mass 
ratio, we have found that the improvement noted by Vec- 
chio holds rather broadly. The degeneracy breaking due 
to precession physics is a rather robust phenomenon. 

Two conclusions from this work are particularly impor- 
tant with regard to the astrophysical reach of future LISA 
measurements. The first is that modeling spin-precession 
physics makes it possible to determine the magnitudes of 
the spins of the black holes which constitute the binary. 
If the spins are rapid, they can be measured quite ac- 
curately (as good as 0.1% accuracy for high spin, low 
redshift systems) due to the strong modulation imposed 
on the signal by their interaction. Coupled with the fact 
that the black hole masses can likewise be measured with 
good precision, this suggests that LISA will be a valuable 
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mi (M©) 


m 2 (M©) 


2a (arcmin) 
(no precession) 


2a (arcmin) 
(precession) 


2b (arcmin) 
(no precession) 


2b (arcmin) 
(precession) 


ADl/Dl 
(no precession) 


ADl/Dl 
(precession) 


10 5 


10 5 


133 


27.3 


84.7 


13.3 


0.0193 


0.00398 


3 x 10 5 


10 s 


115 


16.9 


72.6 


7.33 


0.0165 


0.00240 


3 x 10 5 


3 x 10 5 


101 


23.3 


62.8 


11.8 


0.0143 


0.00357 


10 6 


10 5 


105 


27.2 


65.1 


6.62 


0.0149 


0.00320 


10 6 


3 x 10 5 


93.1 


31.3 


57.5 


13.2 


0.0132 


0.00393 


10 6 


10 6 


90.1 


40.2 


54.1 


21.9 


0.0125 


0.00560 


3 x 10 6 


3 x 10 5 


95.0 


34.1 


57.3 


9.20 


0.0135 


0.00376 


3 x 10 6 


10 6 


102 


32.3 


56.0 


14.7 


0.0135 


0.00419 


3 x 10 6 


3 x 10 6 


135 


43.3 


68.5 


22.3 


0.0182 


0.00689 


10 7 


10 6 


149 


37.6 


75.2 


12.2 


0.0200 


0.00457 


10 7 


3 x 10 6 


238 


42.1 


119 


19.0 


0.0322 


0.00610 


10 7 


10 7 


466 


81.3 


232 


38.6 


0.0636 


0.0250 



TABLE IV: Median errors in extrinsic quantities for 10 4 binaries of various masses at z = 1, including comparisons with the 
"no precession" case. Note that the given major axis and minor axis are the medians for each data set and do not correspond 
to the same binary. However, they still represent an average sky position error ellipse in the following sense: VtyoE, calculated 
using the median values of 2a and 26, differs in most cases by less than 10% from the median value of \/Af2jv calculated from 
the covariance matrix and (|4.13|l (except at more extreme mass ratios — when mxjmi = 10, the difference can be 25%). 



mi (M ) 


m 2 (M©) 


2a (arcmin) 
(no precession) 


2a (arcmin) 
(precession) 


2b (arcmin) 
(no precession) 


2b (arcmin) 
(precession) 


AD L /D L 
(no precession) 


AD L /D L 
(precession) 


10 5 


10 5 


432 


81.0 


271 


40.8 


0.0617 


0.0123 


3 x 10 5 


10 5 


389 


92.5 


242 


39.5 


0.0551 


0.0126 


3 x 10 5 


3 x 10 5 


356 


142 


220 


75.7 


0.0502 


0.0201 


10 6 


10 5 


379 


141 


233 


36.6 


0.0550 


0.0155 


10 6 


3 x 10 5 


359 


129 


215 


56.7 


0.0500 


0.0161 


10 6 


10 6 


416 


158 


224 


84.3 


0.0556 


0.0237 


3 x 10 6 


3 x 10 5 


425 


132 


233 


40.3 


0.0568 


0.0153 


3 x 10 6 


10 6 


599 


142 


302 


64.6 


0.0809 


0.0193 


3 x 10 6 


3 x 10 6 


990 


224 


494 


111 


0.134 


0.0422 


10 7 


10 6 


1320 


206 


648 


78.5 


0.178 


0.0293 


10 7 


3 x 10 6 


2380 


297 


1180 


152 


0.326 


0.0805 


10 7 


10 7 


6820 


2000 


3390 


583 


0.935 


2.41 



TABLE V: Median errors in extrinsic quantities for 10 4 binaries of various masses at z = 3. 



tool for tracking the evolution of both mass and spin 
over cosmic time. Such observations could provide a di- 
rect window into the growth of cosmological structures. 
Measuring spin may also make it possible to indirectly 
test the black hole area theorem [62[. The requirement 
that black hole area can only grow implies a consistency 
relation between the initial and final masses and spins. 
By measuring the initial masses and spins through the 
inspiral, and the mass and spin of the merged remnant 
hole through the ringdown waves [38], [3^ | , we can check 
this consistency relation in a manner analogous to the 
mass loss test proposed in (63[. We intend to investigate 
whether this test is feasible in future work. 



Second, we confirm Vecchio's result that precession 
breaks degeneracies between the angles which determine 
a binary's orientation and its position on the sky, improv- 
ing the accuracy with which sky position can be fixed 
using GWs alone. At low redshift (z ~ 1), we find that 
sources can be localized to within an ellipse whose major 
axis is typically a few x 10 arcminutes across and whose 
minor axis is typically a factor ^2 — 4 smaller. This is 
small enough that searching the GW pixel for an electro- 
magnetic counterpart to the merger event should not be 
too arduous a task [35]. For mergers at higher redshift, 
the waves weaken and the source is not so well localized. 
The field which would need to be searched for sources at 
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mi (M©) 


m 2 (M©) 


2a (arcmin) 
(no precession) 


2a (arcmin) 
(precession) 


2b (arcmin) 
(no precession) 


2b (arcmin) 
(precession) 


ADl/Dl 
(no precession) 


ADl/Dl 
(precession) 


10 5 


10 5 


729 


169 


456 


85.7 


0.104 


0.0260 


3 x 10 5 


10 s 


676 


217 


419 


95.8 


0.0957 


0.0284 


3 x 10 5 


3 x 10 5 


650 


295 


395 


161 


0.0917 


0.0409 


10 6 


10 5 


686 


248 


416 


66.8 


0.0983 


0.0273 


10 6 


3 x 10 5 


716 


233 


404 


101 


0.0961 


0.0294 


10 6 


10 6 


976 


315 


497 


162 


0.132 


0.0501 


3 x 10 6 


3 x 10 5 


986 


265 


507 


86.4 


0.133 


0.0318 


3 x 10 6 


10 6 


1620 


304 


810 


139 


0.220 


0.0436 


3 x 10 6 


3 x 10 6 


2930 


538 


1460 


260 


0.400 


0.140 


10 7 


10 6 


5080 


577 


2480 


290 


0.689 


0.124 


10 7 


3 x 10 6 


10500 


1720 


5130 


621 


1.42 


1.24 


10 7 


10 7 


75500 


180000 


35000 


29600 


10.3 


377 



TABLE VI: Median errors in extrinsic quantities for 10 4 binaries of various masses at z = 5. Again, the results for the highest 
masses are essentially meaningless-the parameters are completely undetermined. 



z ~ 3 — 5 is typically a few degrees across in the long 
axis and tens of arcminutes to a degree or two in the 
short direction — a rather more difficult challenge, but 
not hopeless. We intend to more thoroughly investigate 
the nature of localization with spin precession, including 
how the pixel evolves with observation time up to final 
merger, in future work. 

As mentioned in the Introduction, our analysis makes 
many assumptions and approximations which are likely 
to affect our results; a goal of future work will be to lift 
these approximations. One major concern is the Gaus- 
sian approximation we have taken to the likelihood func- 
tion. As already discussed, this approximation is known 
to be good when the SNR is "large" 8, 28]; however, it is 
not apparent what large really means, particularly given 
that we are fitting for 15 parameters. Lifting this simpli- 
fying approximation can be done by simply computing 
the likelihood function (|3.9p directly and examining how 
well parameters are thereby determined. In the context 
of GW measurements, Markov Chain-Monte Carlo tech- 
niques have been investigated and found to be very useful 
[64 l65l (6(| ; the application of these techniques to LISA 
measurement problems is now being rather actively in- 
vestigated [srl. [gsL [69l] . 

Because we have used the Gaussian approximation 
(among other simplifications taken in this analysis), 
we cannot claim that this analysis gives a definitive 
statement about the accuracy with which LISA could 
measure binary black hole source parameters. However, 
it is certainly indicative of the accuracy which we 



expect LISA to achieve. In particular, we are confident 
that the trends we have seen as parameters are varied 
(e.g., masses, redshift, spin magnitude) are robust. 
Most importantly, it is very clear that the influence of 
spin-induced precession upon the measured waveform 
allows parameters to be measured to greater accuracy 
than before. 
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Erratum: Measuring coalescing massive binary black holes with gravitational waves: 
The impact of spin-induced precession [Phys. Rev. D 74, 122001 (2006)] 



Ryan N. Lang and Scott A. Hughes 



PACS numbers: 04.80.Nn, 04.25.Nx, 04.30.Db, 99.10.Cd 



Since the publication of this paper, we found a rather silly error in the code used for this analysis. We use the 
variable cos 9 to describe an event's location on the sky and to specify its orbital angular momentum vector. In several 
places, we need to convert to sin 9, which we do using the rule sin (9 = \/l — cos 2 9. In two instances (out of eight), we 
failed to square the cost? under the square root. This error impacted both the code that includes black hole precession 
physics and the "no precession" code. 

In fixing this error, we find that our qualitative results are unchanged. However, almost all of our results are 
quantitatively impacted at least slightly. Below we show the corrected versions of Tables I - VI, which give the 
median errors in intrinsic signal parameters (masses and spins; Tables I - III) and extrinsic signal parameters (sky 
position and distance; Tables IV - VI) for a range of masses and redshifts. 
^"v Comparison of these tables to the published versions shows the impact of the coding error on our results. Considering 
the "precession" results first, we find that the corrected errors in mi and 777,2 are typically larger by about 10% — 20% 
than previously. The spin errors follow the same pattern, although in some cases the differences reach ~ 40%. The 
errors to the chirp mass M. are typically worsened by only a few percent (and actually improve a bit for larger masses). 
The reduced mass fi behaves similarly to mi and 7712, showing an accuracy degradation of 10% — 20%. 

Turning next to the "no precession" results, we find that the errors in chirp and reduced mass are typically slightly 
iy) , better than our previous results by a few percent. Taken together with the fact that correcting this bug often degrades 
^ ■ our precession results a bit, we see that the improvement in parameter estimation is not always quite as good as we 
had originally claimed. It is still, however, a significant improvement. For example, precession physics improves the 
accuracy with which the reduced mass can be determined by one to several orders of magnitude. 

Because of the nature of the coding error, errors in the extrinsic parameters were more strongly affected than errors 
in the intrinsic parameters. With the precession code corrected, the major axis of the sky position error ellipse, 2a, 
and errors in the luminosity distance Dl increase by almost a factor of 2 in the worst cases. The "no precession" 
results change by similar factors. So, on average, the improvement due to precession is essentially unaffected by this 
""^ ■ bug for both 2a and Dl- 

The story is somewhat different for the minor axis of the sky position error ellipse, 26. Its accuracy is also degraded, 
but only by a factor ~ 1.5 or so (in the worst cases). In all cases, we find that this bug has less impact on 26 than 
it does on 2a. Our corrected precession code thus indicates that the ratio 2a /2b is more extreme, approaching 4 in 
• • . large mass ratio binaries with enhanced precession effects. Interestingly, we find that 26 for the "no precession" case 
is degraded by as much as a factor ~ 3. Correcting this bug therefore indicates that the minor axis of the localization 
, ellipse is improved by precession even more than we previously reported. 

The changes in the intrinsic parameter accuracy are so slight that the the distributions shown in Figs. 2-6 are 
not noticeably affected. However, the distributions of the extrinsic parameters change significantly enough that we 
present new versions of Figs. 7-9 below. Despite the various changes in shape, the distribution of 26 still contains a 
long tail of small errors when precession effects are included. 

An updated version of our manuscript, including these new results and with small changes in the summary text to 
account for them, has been posted to arxiv.org. 
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TABLE I: Median errors in intrinsic quantities for 10 4 binaries of various masses at z = 1, including comparisons with the "no 
precession" case where possible. We have omitted the errors in chirp mass and reduced mass for equal mass binaries because 
that parameterization of the waveform fails the Gaussian approximation at those points. 
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m 2 (Mq) 
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TABLE II: Median errors in intrinsic quantities for 10 4 binaries of various masses at z = 3. 
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mi (M ) 


m 2 (M Q ) 
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Am,2/m,2 




AX2 


AM/M 
(no precession) 
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0.0156 
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0.00427 


0.000289 


0.469 J 


0.00596 


10 6 


IO 6 


0.0424 


0.0423 


0.197 


0.200 










3 x 10 6 


3 x IO 5 


0.0161 


0.0117 


0.0158 


0.0808 


0.0115 


0.00103 


0.643 


0.00922 


3 x 10 6 


10 6 


0.0576 


0.0475 


0.0606 


0.107 


0.0108 


0.00265 


0.635 


0.0214 


3 x 10 6 


3 x 10 6 
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10 7 


IO 7 


2280 
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TABLE III: Median errors in intrinsic quantities for 10 4 binaries of various masses at z — 5. The results for the highest masses 
are meaningless — the parameters are completely undetermined. 



mi (M©) 


m 2 (M Q ) 


2a (arcmin) 
(no precession) 


2a (arcmin) 
(precession) 


2b (arcmin) 
(no precession) 


2b (arcmin) 
(precession) 


AD L /D L 
(no precession) 


AD L /D L 
(precession) 


10 5 


10 5 


133 


27.3 


84.7 


13.3 


0.0193 


0.00398 


3 x 10 5 


10 5 


115 


16.9 


72.6 


7.33 


0.0165 


0.00240 


3 x 10 5 


3 x 10 5 


101 


23.3 


62.8 


11.8 


0.0143 


0.00357 


10 6 


10 5 


105 


27.2 


65.1 


6.62 


0.0149 


0.00320 


10 6 


3 x 10 5 
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31.3 


57.5 


13.2 


0.0132 


0.00393 
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40.2 


54.1 
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57.3 
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0.00419 
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22.3 
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0.00689 


10 7 


10 6 
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37.6 


75.2 


12.2 
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0.00457 


10 7 


3 x 10 6 
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42.1 


119 


19.0 


0.0322 


0.00610 


10 7 


10 7 


466 


81.3 


232 


38.6 


0.0636 


0.0250 



TABLE IV: Median errors in extrinsic quantities for 10 4 binaries of various masses at z = 1, including comparisons with the 
"no precession" case. Note that the given major axis and minor axis are the medians for each data set and do not correspond 
to the same binary. However, they still represent an average sky position error ellipse in the following sense: V nab, calculated 
using the median values of 2a and 26, differs in most cases by less than 10% from the median value of ^/AQn calculated from 
the covariance matrix and (4.13) (except at more extreme mass ratios — when m\/m2 = 10, the difference can be 25%). 
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mi (M Q ) 


m 2 (M ) 


2a (arcmin) 
(no precession) 


2a (arcmin) 
(precession) 


2b (arcmin) 
(no precession) 


2b (arcmin) 
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10 5 
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81.0 
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40.8 


0.0617 


0.0123 


3 x 10 5 
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92.5 


242 


39.5 


0.0551 


0.0126 


3 x 10 5 


3 x 10 5 


356 


142 


220 


75.7 


0.0502 


0.0201 


10 6 


10 5 


379 


141 


233 


36.6 


0.0550 


0.0155 


10 6 


3 x 10 5 


359 


129 


215 


56.7 


0.0500 


0.0161 


10 6 


10 6 


416 


158 


224 


84.3 


0.0556 


0.0237 


3 x 10 6 


3 x 10 5 


425 


132 


233 


40.3 


0.0568 


0.0153 


3 x 10 6 


10 6 


599 


142 


302 


64.6 


0.0809 


0.0193 


3 x 10 6 


3 x 10 6 


990 


224 


494 


111 


0.134 


0.0422 


10 7 


10 6 


1320 


206 


648 


78.5 


0.178 


0.0293 


10 7 


3 x 10 6 


2380 


297 


1180 


152 


0.326 


0.0805 


10 7 


10 7 


6820 


2000 


3390 


583 


0.935 


2.41 



TABLE V: Median errors in extrinsic quantities for 10 4 binaries of various masses at z — 3. 



mi (M Q ) 


m 2 (M ) 


2a (arcmin) 
(no precession) 


2a (arcmin) 
(precession) 


26 (arcmin) 
(no precession) 


2b (arcmin) 
(precession) 


ADl/Dl 
(no precession) 


ADl/D l 
(precession) 


10 5 


10 5 


729 


169 


456 


85.7 


0.104 


0.0260 


3 x 10 5 


10 5 


676 


217 


419 


95.8 


0.0957 


0.0284 


3 x 10 5 


3 x 10 5 


650 


295 


395 


161 


0.0917 


0.0409 


10 6 


10 5 


686 


248 


416 


66.8 


0.0983 


0.0273 


10 6 


3 x 10 5 


716 


233 


404 


101 


0.0961 


0.0294 


10 6 


10 6 


976 


315 


497 


162 


0.132 


0.0501 


3 x 10 6 


3 x 10 5 


986 


265 


507 


86.4 


0.133 


0.0318 


3 x 10 6 


10 6 


1620 


304 


810 


139 


0.220 


0.0436 


3 x 10 6 


3 x 10 6 


2930 


538 


1460 


260 


0.400 


0.140 


10 7 


10 6 


5080 


577 


2480 


290 


0.689 


0.124 


10 7 


3 x 10 6 


10500 


1720 


5130 


621 


1.42 


1.24 


10 7 


10 7 


75500 


180000 


35000 


29600 


10.3 


377 



TABLE VI: Median errors in extrinsic quantities for 10 4 binaries of various masses at z = 5. Again, the results for the highest 
masses are essentially meaningless-the parameters are completely undetermined. 
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FIG. 7: Distribution of the major axis of the sky position error ellipse, 2a, for 10 4 binaries with mi = 1O 6 M0 and 7712 = 
3 x 1O J M0 at z = 1. The dashed line is the precession-free calculation; the solid line includes precession. Sky position, as an 
extrinsic parameter, is improved somewhat indirectly by precession; therefore, the improvement is less than for the masses. 
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FIG. 8: Distribution of the minor axis of the sky position error ellipse, 2b, for 10 4 binaries with mi = 1O 6 M0 and m2 = 3x 10 5 M© 
at z = 1. The dashed line is the precession-free calculation; the solid line includes precession. 
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FIG. 9: Distribution of errors in the luminosity distance for 10 4 binaries with mi = 1O 6 M and m.2 = 3 x 1O 5 M at z = 1. 
The dashed line is the precession-free calculation; the solid line includes precession. 



Erratum: Measuring coalescing massive binary black holes with gravitational waves: 
The impact of spin-induced precession [Phys. Rev. D 74, 122001 (2006)] 



Ryan N. Lang and Scott A. Hughes 



PACS numbers: 04.80.Nn, 04.25.Nx, 04.30.Db, 99.10.Cd 



Wc have discovered an important typographical error in equations (4.2) and (4.3) of this paper. In each equation, 
a factor of 2 was omitted. The corrected equations are: 




